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We perform a lattice investigation of QCD with three colors and 2 flavors of Dirac (staggered) 
fermions in the adjoint representation, defined on a Ad space with one spatial dimension compactified, 
and study the phase structure of the theory as a function of the size Lc of the compactified dimension. 
We show that four different phases take place, corresponding to different realizations of center 
symmetry: two center symmetric phases, for large or small values of Lc, separated by two phases 
in which center symmetry is broken in two different ways; the dependence of these results on the 
quark mass is discussed. We study also chiral properties and how they are affected by the different 
realizations of center symmetry; chiral symmetry, in particular, stays spontaneously broken at the 
phase transitions and may be restored at much lower values of the compactification radius. Our 
results could be relevant to a recently proposed conjecture Q of volume indepedence of QCD with 
adjoint fermions in the large Nc limit. 
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I. INTRODUCTION 

Extracting the physics of large Nc QCD from simula- 
tions or analytic computations performed on small vol- 
umes, down to a single site, has been the aim of large 
theoretical efforts since the original work by Eguchi and 
Kawai However, the expectation that physics on a 
d-dimensional Li x . . . Ld periodic lattice should be in- 
dependent, as iVc ^ 00, of the lattice sizes Li, . . . Ld, is 
justified if no transition occurs as the sizes Li, . . . Ld are 
changed: this is not true in ordinary Yang-Mills theo- 
ries in d > 2, which undergo phase transitions as one or 
more dimensions are compactified below a given thresh- 
old, corresponding to the spontaneous breaking of the 
center symmetry in that direction. One or more Wilson 
(Polyakov) lines in the compactified directions get a non- 
zero expectation value, i.e. one or more deconfining-like 
transitions take place (see Ref. Q for a review on the 
subject). 

Large Nc QCD with Nf fermions in the antisymmetric 
representation, or QCD(AS), has attracted lot of theoret- 
ical interest in the last few years, as a different Nc — > 00 
limit of ordinary QCD: it is indeed equivalent to ordinary 
QCD for Nc — 3. Orientifold planar equivalence states 
that QCD(AS) is equivalent, in the large Nc limit and 
in the charge-even sector of the theory, to QCD with 
fermions in the adjoint representation, QCD(Adj) [1]. 
Such equivalence is guaranteed if the charge conjugation 
symmetry is not spontaneously broken in QCD (AS) 4]: 
that is true for sufficiently large volumes, but fails in 
presence of a compactified dimension, below a typical 



compactification radius [B!, ^ . 

Recently it has been made the hypothesis that in- 
stead QCD with one or more adjoint fermions would 
persist in the confined phase (i.e. no phase transition 
would occur) as one of the dimensions is compactified, 
if the boundary conditions for fermions are periodic in 
that direction (i.e. if the compactified dimension is not 
temperature-like). That would imply physics be volume 
independent in the large Nc limit, and thus would open 
the possibility of studying QCD(Adj) down to arbitrarily 
small volumes, finally connecting information, via orien- 
tifold planar equivalence, to the physical properties of 
QCD(AS) on large volumes. That would be a break- 
through in the study of large Nc QCD, with possible im- 
portant implications also for the knowledge of real QCD. 
ft is therefore of great importance to work on such con- 
jecture, providing evidence from numerical lattice simu- 
lations. 

The argument given in Ref. [7] is based on the analysis 
of the 1-loop effective potential for the Polyakov (Wilson) 
line, n, taken along a compactified dimension of length 
Lc, in presence of Nf {Nf/2) flavors of massless Majo- 
rana (Dirac) adjoint fermions, which is given, up to a 
constant, by 
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where the —/+ sign occurs for periodic/antiperiodic 
boundary conditions (p.b.c./a.b.c.) in the given di- 
rection, 6**^^ , j — l,...Nc are the eigenvalues of the 
Polyakov line, which satisfy the constraint Ylj e"^^ = 1, 
and [x] = {x mod 27r). 

At arbitrarily small compactified dimension, the 1-loop 
potential gives all relevant information. For Nf = or for 
a.b.c, the potential implies attraction among the eigen- 
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values: indeed it is minimized when they coincide, mean- 
ing that the Polyakov line alignes itself along some center 
element and center symmetry is spontaneously broken. 
Instead, for Nf > 1 and p.b.c. the potential changes sign 
and becomes repulsive: the Polyakov line is disordered 
in this case and center symmetry is not broken, exactly 
as in the infinite volume, confined phase. The potential 
in Eq. ([T]) is given for zero quark mass m and correc- 
tions are expected for finite m, however results should 
not change as finite quark masses are switched on, and 
indeed it has been shown in Ref. @ that the disordered 
confined phase is always preferred for small enough (com- 
pared to Aqcd) quark masses. Analogous results are ob- 
tained when studying the I-loop potential in the lattice 
regularized theory II], while contrasting results (i.e. a 
center symmetry breaking I-loop potential) obtained in 
a three-dimensional reduced model jl^ should be inter- 
preted with care and may be due to the absence of cer- 
tain center stabilizing relevant operators which should in 
principle be included in the reduced model [ij]. Finally, 
a similar effect is obtained in the pure SU(A^c) gauge 
theory by adding a trace deformation (i.e. proportional 
to the adjoint trace of the Polyakov line) to the usual 
Yang-Mills action 

An essential requirement for volume independence in 
the large Nc limit is that nothing, i.e. no phase transi- 
tion, happens as the compactified dimension spans from 
infinity to the arbitrarily small, perturbative regime, i.e. 
that the system stays always confined and does not ex- 
plore any other phase in between. That can be verified 
by numerical lattice simulations. While numerical results 
exist for pure SU(Afc) gauge theories with trace defor- 
mations flOl, no result still exists concerning QCD with 
dynamical adjoint fermions. 

In the present paper we consider the example of Nc = 3 
with 2 flavors of Dirac adjoint fermions, providing nu- 
merical evidence that the requirement above is highly 
non-trivial. The system that we have studied undergoes 
several phase transitions when shrinking from the large 
volume to the small volume confined regimes. Our study 
has been performed for finite values of the quark masses: 
the fact that these phase transitions may disappear in the 
chiral limit is not excluded, but our present data show 
that this is non-trivial. 

While our findings could be relevant just to the par- 
ticular case studied here and to the discretization setup 
used, they are a warning claiming for further studies with 
different values of Nc and approaching the continuum 
limit. 

The paper is organized as follows. In Section II we 
present details about lattice QCD with adjoint fermions 
and about our numerical simulations. In Section III we 
present our numerical results, concerning the realization 
of both center symmetry and chiral symmetry as a func- 
tion of the size of the compactified dimension. In Section 
IV we discuss our conclusions. 



II. QCD (AD J) ON THE LATTICE AND 
SIMULATION DETAILS 

Fermions in the adjoint representation of SU{3) have 
8 color degrees of freedom and can be described by 3 x 3 
hermitian traceless matrices: 

Q = Q^Xa (2) 

where Aa are the Gell-Mann's matrixes. The elemen- 
tary parallel trasports acting on them are given by the 
8— dimensional t/jg) representation (which is real) of the 
gauge links 

C^(8).,. = ^Tr(A'^C/(3),,A''[/;3) J (3) 

where z is a lattice site. The full discretized lattice action 
used in our investigation is given by 

S^Sg [C/(3)] + E M t (4) 

Sg is the standard pure gauge plaquette action with links 
in the 3-dimensional representation 

^a = /^E(l-]^Trno) 

□ \ c / 

where (3 is the inverse gauge coupling (3 = 2Nc/gQ and 
the sum is over all elementary plaquette operators IIq of 
the lattice. M is the fermionic matrix, which has been 
chosen according to the standard staggered fermion dis- 
cretization: 

^ 4 

iy=l 

Notice that, since gauge links acting on adjoint fermions 
are real matrixes, M°^^ is real as well. The functional 
integral, describing 2 flavors of Dirac adjoint fermions in 
the continuum limit, is then given by 

Z = y"x>t/e"^^(detM[C/(8)])i/2 (g) 

Due to the reality of ill, that can be expressed in terms 
of a real pseudofcrmion field living only on even sites of 
the lattice as follows 

Z = J VUV^ce-^" exp (M*M)-i4>e) (7) 

In this form the functional integral distribution can be 
easily sampled by a standard Hybrid Monte Carlo algo- 
rithm. In particular we have adopted the so-called <& 
algorithm [13l |. 

We shall consider in the following a system in which 
one dimension is compactified and has size Lc- Adjoint 
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parallel transports are blind to the center of the gauge 
group, hence the presence of adjoint fermions does not 
break explicitely center symmetry in the compactified di- 
rection. The fundamental Polyakov line L(3) taken along 
the compactified dimension is not invariant under the 
same symmetry, and is therefore an exact order parame- 
ter for its spontaneous breaking. 

If fermions are given a.b.c. in the compactified dimen- 
sion, the Feynman path integral can be interpreted as 
the partition function at temperature T = 1/Lc- The 
spontaneous breaking of center symmetry is then associ- 
ated to deconfinement. Various numerical studies of fi- 
nite temperature QCD(Adj) have been performed in the 
past, mostly aimed at studying the relation between de- 
confinement and chiral symmetry restoration [T^j fisl . [l6| . 

In our investigation we have considered fermions with 
p.b.c. along the compactified dimension. This the- 
ory cannot be interpreted as a thermodynamical system 
(apart from the case of infinite fermion mass), but must 
considered simply as a system with a shorter, compacti- 
fied spatial dimension. Nevertheless we shall refer to the 
phase where the fundamental Polyakov line acquires a 
non-zero expectation value, oriented along an element of 
the center of the gauge group, as a "deconfined" phase. 

In most of our numerical simulations we have kept a 
fixed number of lattice sites along the compactified di- 
rection, and varied its physical size by tuning the lattice 
spacing a via the inverse gauge coupling [3. Asymptotic 
freedom implies that a — > as /3 — > cxd, hence the com- 
pactified dimension is made shorter and shorter as (3 is in- 
creased. Since we are mostly interested in studying which 
phases are explored by the system as the size of the com- 
pactified dimension is changed from large to small values, 
and not in the exact location of the possible phase tran- 
sitions, it is enough for our purposes to look at the phase 
structure of theory as a function of [3 at fixed number of 
lattice sites. However we shall try to get an estimate of 
physical scale ratios, by means of the 2-loop perturba- 
tive /3-function, describing the asymptotic dependence of 
the lattice spacing a on (3, or alternatively by performing 
numerical simulations also at a fixed value of the lattice 
spacing, in which is tuned by changing the number of 
lattice sites in the compactified dimension. 

It is essential to repeat our investigation for different 
values of the quark mass. An infinite quark mass corre- 
sponds to the pure gauge theory, where it is well known 
that two different phases, deconfined or confined, take 
place respectively for short or large values of the com- 
pactified dimension; such phases are separated by a first 
order transition 17]. For small enough quark masses one 
expects instead center symmetric phases both for a very 
large or a very short size of the compactified dimension, 
with or without different phases in between. It is then 
clear that a highly non-trivial phase structure must be 
found in the quark mass - plane. For this reason 
we have performed numerical simulations at five differ- 
ent values of the bare quark mass, am — 0.50, 0.10, 0.05, 
0.02 and 0.01. 



In all cases we have chosen molecular dynamics trajec- 
tories of length 0.5, with typical integration steps 5t rang- 
ing from 0.025 (for am = 0.50) to 0.005 (for am = 0.01), 
so as to keep an acceptance rate around 80%. For each 
parameter choice (am, (3, L^) we have collected between 
2k and 20k trajectories. Numerical simulations have been 
performed on the APEmille machine in Pisa and the 
apeNEXT facility in Rome. 

III. NUMERICAL RESULTS 

We start by presenting results obtained on a lattice 
16^^ X Lc, with Lj, = 4 and p.b.c. for all fields. The 
shorter direction plays the role of the compactified di- 
mension of length Lc = Lca{f3,m). In order to distin- 
guish the various phases corresponding to different pos- 
sible realizations of center symmetry, we have studied 
the expectation value and distribution of the trace of the 
Polyakov line, L(3) , taken along the shorter direction and 
averaged over the orthogonal 3d space, as a function of 
/3 and am. 

In Figs. [1] and [5] we show scatter plots of the distri- 
bution of L(3') in the complex plane, for am = 0.1 and 
am — 0.02 respectively, and for different values of f3. In 
each case four different phases are clearly distinguishable: 

• In a first phase at low /3 values, correponding to 
large values of Lc, i(3) is distributed around zero 
with a vanishing expectation value. That corre- 
sponds to the usual confined, center symmetric 
phase of QCD(Adj). 

9 As (3 increases (Lc is made shorter), the system 
passes into a phase with a broken center symmetry, 
with L(3) aligned along one of the center elements 
of SU(3). We call this phase deconfined, since it is 
analogous to the usual high temperature phase of 
QCD. 

• At very high values of f3 (very small Lc) the system 
is again in a center symmetric, confined-like phase, 
in agreement with the prediction of the one-loop 
potential in Eq. Since this phase seems to be 
distinguished and separated from the low temper- 
ature confined phase, we refer to it as re-confined. 

• Finally, a new phase is present for values of /3 (Lc) 
between the deconfined and the re-confined phase, 
in which center symmetry is broken by a non-zero 
expectation value of L(3) , which is however ori- 
ented along directions rotated roughly by an an- 
gle TT, in the complex plane, with respect to those 
corresponding to center elements. This phase is 
completely analogous to that found by the authors 
of Ref. [lol in their study of pure gauge SU(7Vc) 
gauge theories with adjoint Polyakov line deforma- 
tions, and predicted again by the same authors, on 
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the basis of the 1-loop effective potential, also in 
the case of QCD(Adj) with p.b.c. and for some 
range of m and Lc [1] . We shall refer to this phase 
as a split or skewed phase, following the convention 
in Ref. M. 
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FIG. 1: Scatter plots of the distribution of the trace of the 
Polyakov line, 1/(3), in the complex plane, for various values 
of f3 and at a fixed quark mass am = 0.10 on a 16"^ x 4 lattice. 



new center symmetric phase with a small average mod- 
ulus for i(3)- In the case of am = 0.5 only two phases 
(confined and deconfined) are visible in Fig. [31 simply be- 
cause of the limited range of P shown in the figure. The 
evidence for the existence of a split (skewed) phase is 
weaker at the smallest quark mass explored, am = 0.01. 
Typical time histories for the real and imaginary parts 
of L(3) are shown in Fig. 01 in particular for the confined 
and deconfined phase at the lowest quark mass explored, 
am = 0.01. 
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FIG. 2: Same as in Fig. [l] but for am = 0.02. 

The transitions among the different four phases are 
also clear from Fig. [31 where we report the average value 
of |L(3)| as a function of (3 for the different quark masses 
explored. In most cases the average modulus is small at 
low /3 (confined) phase; then starts growing at a critical 
value Pc/d separating the confined from the deconfined 
phase; a second critical value Pd/s is reached, separating 
the deconfined from the split (skewed) phase, at which 
the modulus drops roughly by a factor 3; finally a third 
critical value Ps/r is met where the system goes into the 



FIG. 3: Average value of the Polyakov line modulus, on a 
16^ X 4 lattice, as a function of f3 and for different bare quark 
masses. 




FIG. 4: Time histories (in units of Molecular Dynamics tra- 
jectories and including part of the thermalization) of the real 
and imaginary part of 1/(3), for am = 0.01 and two different 
P values in the confined and deconfined phase respectively. 

It is interesting to notice that in each different phase 
the Polyakov line modulus has a very small dependence 
on the quark mass, much weaker than what observed in 
the case of finite temperature QCD(Adj) (see for instance 
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TABLE I: Critical values of /3 as a function of the bare quark 
mass for the transition from the confined to the deconfined 
phase {/3c/d), from the deconfined to the split (skewed) phase 
(/3(j/s), and from the split to the weak coupling confined (re- 



confined) phase (/3, 



s/c) 



In the case of am = 0.01 the split 



phase is absent and (5^/^ actually indicates the transition from 
the deconfined to the re-confined phase. 
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FIG. 5: Adjoint Polyakov loop compared to the modulus of 
the fundamental loop (divided by a factor ten to fit in the 
figure), as a function of /3 for am — 0.05 on a 16"^ x 4 lattice. 



results reported in Ref. [Tj, [1^). In the present case the 
main effect of changing the quark mass seems that of 
moving the location of the phase transitions. 

Locations of the critical couplings as a function of the 
bare quark mass are reported in Table IH Regarding 
the order of the phase transitions, we have indications, 
coming from abrupt jumps of observables or presence of 
mctastabilities, that they are all first order, at least for 
am > 0.05. For the smallest quark masses, am = 0.01 
and am = 0.02, transitions look smoother. A deeper 
investigation and a finite size scaling analysis would be 
needed to reach a definite conclusion: that goes beyond 
our present purposes. However we notice that true phase 
transitions (continuous or discontinuous) should be found 
at points where the realization of the (exact) center sym- 
metry changes. 

Let us now discuss our results. The confining be- 
haviour observed at large Lc (low P) is in agreement 
with the usual infinite volume behaviour of QCD(Adj). 
The confining behaviour observed at small Lc (high /3) 
is in agreement with the prediction of the one-loop ef- 



fective potential in Eq. fermions contribute with a 
repulsive, disordering interaction for the Polyakov line 
eigenvalues, which is opposite and wins over the ordering 
gluon contribution. What is non-trivial is the behaviour 
observed at intermediate values of where it is clear 
that higher order and finite quark mass corrections to the 
one-loop effective potential play a significant role. In the 
deconfined phase those corrections result in an overall at- 
tractive interaction for the eigenvalues which orders the 
Polyakov line along a center element. In the case of the 
split (skewed) phase instead, the free energy minimum is 
apparently reached as two eigenvalues coincide and the 
third one gets a phase tt with respect to them, as in the 
following SU(3) matrix 




(8) 



That is not a minimum for the potential in Eq. ([T]), 
which is instead minimized (in the case of p.b.c. and for 
Nf > 1) by the traceless matrix dia g(l, e^''^/'^, e~^'^/^), 
however it has been verified in Ref. that mass cor- 
rections to the 1-loop potential may move the minimum 
to group elements like that in Eq. ([5]). The trace of such 
group elements is 1/3 with respect to that of center ele- 
ments: that roughly explains the drop observed for the 
modulus of L(3) at Pd/s- 

It is also interesting to look at the behaviour of the 
adjoint Polyakov loop, which is shown in Fig. [5] for the 
case am = 0.05: it is positive only in the deconfined phase 
and negative outside. Since Tr(8) = |Tr(3) p — 1, the large 
negative values obtained for the average adjoint Polyakov 
loop in the re-confined phase can be interpreted as the 
fundamental loop becoming traceless point by point in 
the weak-coupling limit, in agreement with the 1-loop 
potential in Eq. ([1]), while in the strong coupling regime 
the fundamental trace averages to zero mainly because 
of long range disorder. 

In the case of bare quark mass am — 0.10 we have re- 
peated our runs on a 16'^ x 6 lattice, i.e. taking Lc — 6 
sites in the compactified direction. Results for the mod- 
ulus of L(-3) are reported in Fig. [51 where they are com- 
pared with results obtained on the 16'^ x 4 lattice. The 
phase structure does not change and all the critical cou- 
plings move to larger values: that is consistent with the 
possible existence of a continuum limit for the critical 
values of the compactified dimension Lc = Lca{(3, am). 
In particular the deconfinement critical coupling moves 
hom (3c/d{Lc = 4) = 5.41(1) to Pc/diL = 6) = 5.50(1): 
that is consistent with a constant critical length Lc/d — 
Lca{(3c/d), since using the 2- loop /3-function 
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we get 



a{P,^d{L, = 4)) 
a{Pc/d{Lc = 6) 



1.6(2) 



Pd/s a-nd Ps/r show a larger change when going from 
Lc = 4: to Lc — Q, which cannot be interpreted in terms 
of the 2-loop /3-fmiction. This is due to the fact that we 
are not working at fixed physical quark mass, but instead 
at fixed am, meaning larger and larger m as /3 increases: 
as we increase we have to increase (3 to make a smaller 
and keep Lc = aLc fixed, but at the same time, if am 
is fixed, also m = {am) /a increases meaning that also 
the physical value of Lc changes. While the critical L^/d 
corresponding to deconfinement changes smoothly and 
stays finite as the quenched limit m — > oo is approached, 
so that this effect may be small, the other two transitions 
must disappear in the same limit, meaning a steepest de- 
pendence of the physical critical dimension on m (hence 
on f}). We notice that for Lc — Q the transition from the 
deconfined to the split (skewed) phase is accompanied 
by strong metastabilities, testified by different simula- 
tions at equal values of [3 staying in different phases for 
thousands of trajectories (see Fig. [S] where the results of 
those simulations are reported for (3 = 6.2 and (3 = 6.4). 

In order to get an idea of the scale separation between 
the different phases we have also performed simulations 
at fixed quark mass and ultraviolet (UV) cutoff, and vari- 
able Lc (only even values of Lc can be explored because 
of the staggered fermion discretization). In particular in 
Fig. [7| we show the distribution of L(3) in the complex 

plane and its average modulus for different value of Lc 
on a 16'^ X Lc lattice, at fixed bare quark mass am = 0.1 
and inverse gauge coupling /5 = 5.75. The system is in 
the confined, center symmetric phase for L — 16, where 
(L(3)) = (0.0001(2), 0.0002(2)); it is stiU marginally con- 
fined for Lc = 14, where (i(3)) = (0.0005(2), 0.0001(2)), 
while as the compactified dimension is squeezed further 
a finite Polyakov line expectation value develops already 
for Lc = 12, where (L(3)) = (0.0066(8), -0.0008(8)): this 
is clearly different from zero, even if quite small due to 
strong UV noise. The system stays in the deconfined 
phase till Lc = 6, while it is clearly in the split phase 
for Lc = 4. Finally, for the shortest compactified di- 
mension explored, Lc = 2, it is in the center symmetric, 
re-confined phase. We conclude that, in this case, the 
two center symmetric, confined phases, taking place at 
large or small values of the compactified dimension, are 
not connected: they have boundaries which are separated 
by a scale factor ~ 3-4. 



A. Approaching the chiral limit 

Our results do not support the conjecture about vol- 
ume independence of QCD(Adj), since for all explored 
masses the two center symmetric phases taking place for 
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FIG. 6: Comparison of results obtained for the average value 
of the Polyakov line modulus as a function of /3 at am = 0.10 
and on two different lattices, 16^ x 4 and 16'' x 6. 





1 ' 1 

i L =2 

1 c 


0.1 


-0.1 




1 ' 

i L =4 

1 ^ 


0.1 

-0.1 




1 

L 


= 6 

c — 




i , 1 








1 




-0.1 


0.1 




-0.1 


0.1 




-0.1 





0.1 




1 ' 1 

i L =8_ 

1 c 


0.1 




1 ' 

1 L^,= 10- 


0.1 




1 

L 

1 


= 16_ 




-|-|»-- 







--- 1* 























i , 1 


-0.1 




i , 


-0.1 




i 





FIG. 7: Scatter plots of the distribution of 1/(3) in the com- 
plex plane on a 16^ x Lc lattice at fixed quark mass and UV 
cutoff (/3 — 5.75 at am = 0.10) for various values of the com- 
pactified dimension Lc- 



large and small values of Lc do not correspond to a unique 
confined phase, but are instead separated by different 
phases iir which center symmetry is spontaneously bro- 
ken. However a very important issue is whether this is 
true for arbitrarily small masses. 

While our computational resources do not allow us to 
perform extensive simulations for smaller quark masses, 
we can try to extrapolate to the chiral limit our results for 
the critical values of (3 obtained on the 16^^ x 4 lattice and 
reported in Table HI We have therefore fitted the critical 
values of (3 obtained for each different phase transition 
according to 



Pcrit = f3c„t + ciam + C2{amf 



(11) 
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where indicates the critical coupHng in the chiral 

hmit. 

Results are reported in Fig. [H The split phase could 
disappear in the chiral limit, and already at am = 0.01 
the evidence for its existence is weaker. The values of the 
critical couplings extrapolated to the chiral limit are the 
following: = 5.310(15) for the transition from the 

confined to the deconfined phase, and = 5.444(15) 
for the transition from the deconfined to the split (or 
re-confined) phase. 

These results provide partial evidence for a deconfined 
phase surviving in the chiral limit. The two confined 
phase are not connected but still separated, even in the 
chiral limit, by a scale factor which is estimated, on the 
basis of the 2-loop /3-function given in Eq. ([9]), to be 
roughly 2. We would like to stress that this is only the 
result of an extrapolation to the chiral limit of the critical 
couplings /3c/d and 13,1/^ delimiting the deconfined phase 
in the phase diagram in Fig. [8l there is still no definite 
evidence that the deconfined phase actually extrapolates 
to am = 0. 

In order to get more information on this issue we have 
performed a numerical simulation at a smaller value of 
the quark mass, am = 0.005, and for a single value of 
(3 chosen well within the deconfined region in Fig. [8l 
[3 — 5.40. In Fig. [H] we report the Monte-Carlo time 
histories obtained for the real and imaginary part of the 
Polyakov loop respectively at am — 0.01 and am — 0.005. 
It is apparent from the figure that while for am — 0.01 
the system stays quite stably in the deconfined phase, for 
the lower mass it spends part of the simulation time into a 
phase with a lower value of the Polyakov loop (seemingly 
a split phase). More extensive (not affordable) simula- 
tions would be necessary to clarify the issue, but we can 
at least say that some metastable behaviour is present at 
this mass which may suggest an even less trivial phase 
structure close to the chiral limit. 



B. Chiral properties 
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FIG. 8; Possible sketch of the phase diagram in the f3-am 
plane and approaching the chiral limit, obtained by means 
of quadratic interpolations of the critical couplings as in 
Eq. ((ni- 



am = 0.01 



— Real part 

— Imaginary part 




am = 0.005 



FIG. 9: Time histories, in units of Molecular Dynamics tra- 
jectories, of the real and imaginary part of 1/(3), for f3 = 5.4 
and two different values of am. 



It is interesting to study the fate of chiral symmetry 
and its connection with the different observed realiza- 
tions of center symmetry in this theory. The study of 
finite temperature QCD(Adj) has shown that, unlike the 
case of ordinary QCD, phases may exist in which decon- 
finement is not accompanied by chiral symmetry restora- 
tion [l3,[ll,[l^. In the present case, since a chirally sym- 
metric phase is expected anyway in the weak coupling 
regime (i.e. for short enough Lc), a new exotic phase 
could exist in which both chiral symmetry and center 
symmetry are not spontaneously broken [l9l . [20l | . 

In Fig. [10] we show the behaviour of the chiral con- 
densate as a function of P for different quark masses, as 
obtained on the 16^ x 4 lattice, while in Fig. [TT]we show 
the extrapolations of the same chiral condensate values 
to the chiral limit for different values of f3 which cover all 
the possible different phases described by the Polyakov 



line. Also in this case we have used quadratic extrapola- 
tions 

{iptp) {am) = {'ipip){0) + aiam + 02 (am)2 (12) 

which are shown as dashed lines in Fig. [Til together with 
similar extrapolations, including also a term proportional 
to \J am (continuous lines), analogous to those used in 
Ref. [l|. 

Fig. [Tl] shows that, independently of the extrapola- 
tion used, the chiral condensate always extrapolates to a 
non-zero value, as am — > 0, in all the range of /3 values 
going at least up to /3 = 6.0. That means that the dif- 
ferent transitions, corresponding to different realization 
of center symmetry, do not affect chiral symmetry, which 
remains spontaneously broken. 



8 



I.. V. _ >^ 



5.2 5.4 5.6 5.8 



• m = 0.01 
■ m = 0.02 

* m = 0.05 
T m = 0.10 



6.2 6.4 



6.6 



0.02 



0.01 



-0.01 



-0.02 - 



-0.03 



-0.04 



G-O Adjoint loop 

▼ Condensate difference (p.b.c. - a.b.c.) 




5.4 



5.6 



5.8 6 
P 



6.2 



6.4 



FIG. 10: Chiral condensate as a function of (3 for different 
values of the bare quark mass on the 16'^ x 4 lattice. 



FIG. 12: Comparison between {{ipil))m — {ipij^ym ) (see text) 
and the adjoint Polyakov loop for am = 0.05 on a 16^ X 4 
lattice. 




am 



FIG. 11: Chiral condensate as a function of the quark mass 
and quadratic extrapolations to the chiral limit including 
(continuous) or not including (dashed) a term proportional 
to y/am, for different values of the gauge coupling on the 
16^ X 4 lattice. 



Instead at the highest value explored, (3 = 6.50, there 
is some evidence that the chiral condensate may extrapo- 
late to zero and chiral symmetry be restored: that would 
fit the expectation of a chirally symmetric weak coupling 
phase. However we cannot draw a definite conclusion 
from present results: at this value of (3 the lattice spac- 
ing is very small and even am = 0.01 corresponds to a 
large quark mass, so that an extrapolation to zero quark 
mass may be debatable; moreover one should also prop- 
erly take into account the running of the adimensional 
quark condensate measured on the lattice as the contin- 
uum limit a — > is approached. 

Therefore, the only reasonable conclusion that we can 



draw from our results is that the expected weak coupling 
restoration of chiral symmetry may happen at values of 
the compactified dimension much shorter (at least one 
order of magnitude, as roughly estimated on the basis of 
the 2-loop /3-function) than those at which center sym- 
metry transitions take place. 

It is interesting to ask if the information about the 
change in the realization of center symmetry is tran- 
scripted in some other way into the chiral properties of 
the theory. A possible candidate to give an answer is the 
dual chiral condensate introduced in Ref. [2l[ and defined 
as 

E, = -^ r d^e-^'^{i,^)(^^ (13) 
27r Jo 

where {'ipip)m'^ is the chiral condensate measured with an 
assigned phase exp(z(p) for the fermionic boundary condi- 
tions in the compactified dimension (ip ^ for p.b.c. and 
if = TT for a.b.c), and Si is nothing but its Fourier trans- 
form with respect to the phase. It can be easily shown, 
e.g. by a loop expansion of the fermionic determinant, 
that in the infinite quark mass limit the dual conden- 
sate becomes the usual Polyakov loop (adjoint Polyakov 
loop in this case), hence it is a natural chiral quantity 
which could be sensible also to the realization of cen- 
ter symmetry, as the adjoint loop is (see Fig. O. As a 
very rough estimate of Ei, we have considered the quan- 
tity (('0V')m'' — (■0'0)m''), i.e. the difference in the con- 
densate measured using p.b.c. and a.b.c. respectively 
(that means only a change in the observable, which is 
measured on the same configurations produced with dy- 
namical fermions having p.b.c). In Fig. [12] we report 
results obtained for am = 0.05: as can be appreciated 
this difference follows quite closely the behaviour of the 
adjoint Polyakov loop, changing its sign in the deconfined 
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phase. This is a confirmation, in a theory with a highly 
non-trivial phase structure, that the realization of center 
symmetry is reflected in the way the chiral condensate 
depends on the fermion boundary conditions f23| . 

IV. CONCLUSIONS 

We have studied QCD with three colors, 2 flavors of 
Dirac (staggered) fermions in the adjoint representation 
and one of the spatial dimensions compactified with pe- 
riodic boundary conditions for fermions. We have shown 
that 4 different phases are explored as the size Lc of 
the compactified dimension is varied, corresponding to 
different realizations of center symmetry. In particular 
two center symmetric, confined phases exist, for small 
and large values of ic, separated by two phases (decon- 
fined and split phase) in which center symmetry is spon- 
taneously broken in two different ways. We have some 
partial evidence that the deconfined phase could persist 
and separate the two confined phases also in the chiral 
limit. This conclusion is based on an extrapolation to 
the chiral limit of the critical lines determined for quark 
masses down to am — 0.01. On the other hand a direct 
simulation at a lower mass, am — 0.005, and at /3 = 5.4, 
which is inside the deconfined phase according to the ex- 
trapolations above, has shown some sign of metastabil- 
ity which could signal a less trivial phase structure close 
to the chiral limit. More extensive simulations at lower 
masses, which are not affordable by our present compu- 
tational resources, would be needed to better clarify the 
issue. 

We have investigated the chiral properties of the the- 



ory and have provided evidence that chiral symmetry is 
not restored until very small values of Lc, at least one 
order of magnitude shorter than those at which center 
symmetry transitions take place. On the other hand, we 
have verified that the realization of center symmetry is 
reflected in the way the chiral condensate depends on the 
fermion boundary conditions, thus confirming similar re- 
sults obtained in Ref. (2^ . 

Our results may be relevant for a recently proposed 
conjecture about volume independence of QCD(Adj) in 
the large Nc limit [tI]: in particular the presence of a 
deconfined phase separating the large and small volume 
confined phases would not support such conjecture. Of 
course, apart from the fact that the extrapolation to the 
chiral limit of our findings should be further investigated, 
our results could be relevant just to the particular case 
studied here and to the discretization setup used. It is 
therefore essential to perform further studies. In particu- 
lar one should verify our results as the continuum limit is 
approached, also making use of improved actions or dif- 
ferent discretization setups, based for instance on Wilson 
instead of staggered fermions. Finally, one should con- 
sider that of course the situation could be different for 
larger values of Nc- 
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